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ABSTRACT 


Crack problems can be solved by using the Distributed Dislocation Technique , where under a 
given loading conditions a crack is modeled as an airay of dislocations since the crack faces 
remain traction free, hence making use of this condition an integral equation is developed which 
is singular in nature, solving this integral equation with the help of quadrature rule we compute 
the dislocation density at the ends of crack tip, which in turn is used to compute the stress 
intensity factor at the tip of a crack. 

After the above analysis this technique is applied to the case of Bi-penodic cracks m an infinite 
body. Results of this analysis are presented with the help of graphs. 

Problem of a inclusion-interface-matrix is also solved by using the same method, here a circular 
inclusion, of a given radius is considered, which is subjected to uniform pressure. 

Results obtained by distributed dislocation technique are also compared with the analytical 
results 
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Chapter 1 
Introduction 


1.1 Introduction 

Finite Element Method (FEM) is one of the numerical methods to obtain an approximate solution 
of the engineering problems. Applications include structures, solid mechanics, heat conduction, 
vibrations, elasto-plastic analysis, fluid flow etc. Finite element method is widely used in fracture 
mechanics applications. In FEM the domain of the problem is discretized into a number of 
subdomains, called finite elements, connected with one another at points known as nodes The 
variables of the problem, such as displacements, temperature, velocity, etc. , aie approximated 
piecewise, so that they are represented in each element by simple polynomials. The coefficients 
of the polynomials are determined such that the governing equation and boundary conditions are 
satisfied in a best possible manner. The solution becomes more and more accurate by choosing 
smaller elements and more number of nodal variables. The necessary and sufficient conditions on 
the nodal variables and the shape functions to assure the convergence are well laid out in FEM 
theory. Some of the typical two dimensional elements used to discretize the domain are shown 
in the figure. 

Two-dimensional problems of fracture mechanics are solved through FEM, in general, by using 
the triangular, quadrilateral or isoparametric elements. However there are certain limitations of 
the simple elements to represent large stress gradients. The stress gradient near the crack tip 
is very high, and theoretically, stress becomes singular at the crack tip. To represent such large 
stresses and stress gradients reasonably well, we need to employ very fine mesh near the crack 
tip. This involves large amount of computational complexity. 

Computed value of SIF has an error of 11% with coarse mesh and 5% when a very fine mesh is 
used. 

There are several indirect methods which do take care of high stress gradient near the crack 
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3-node triangle 


4-node quadrilateral 




Figure 1.1. Typical two- dimensional elements 



Figure 1.2: A rectangular element mesh and path of integration for J-Integral 


tip. They make use of some knowledge of crack behaviour in finite element analysis which will 
markedly improve the accuracy of the results The most popular method is the J- Integral 
method. In this method a path is choosen from one crack surface to the other crack surface thus 



surrounding the crack tip. This integral value can be evaluated on a path which is slightly away 
from the crack tip. The resulting stress intensity factor is quite accurate having an error of 3.5%. 
Another method is by employing the energy release rate technique. This involves the evaluation 
of energy m the cracked body for two slightly different configurations. The energy release rate 
G can be written in terms of SIF. This computed SIF has an erroi of 3% 




Quadrilateral element 

Figure 1.3: Quarter point elements around a crack tip 


Another simulation of singularity is by the use of quarter point element technique. In this 
technique usual isoparametric 6-noded triangle or 8-noded isoparametric quadrilateral elements 
are employed. Midside nodes on the two adjacent side are shifted towards the corner node to the 
quarter point locations. For these locations of the mid nodes the stresses becomes singular .In the 
figure observe that all the midside nodes adjacent to the crack tip are at quarter point location. It 
has been demonstrated that SIF found by this method is within 2% of the theoretical solutions. 

The technique presented in this thesis is based on the pioneering work of Eshelby in the late 
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1950’s, and developed by Erdogan, Keer and Mura. The basic idea of this technique is to use 
the superposition of the stress field present in the unfiawed body, together with an unknown 
distribution of ‘strain nuclei’ (Dislocations), chosen so that the crack faces become traction-free. 
The solution used for the stress field for the nucleus is chosen so that other boundary conditions 
are satisfied. The technique is therefore efficient, and may be used to model the evolution of a 
developing crack in two or three dimensions. 

The basic strategy used is to cancel the unwanted tractions appearing along the faces of flaws by 
distributing various kinds of strain nuclei(Dislocations) along the length of flaw in two dimension. 
Dislocations are the point source of strain (or stress). The problem is solved in three parts, first 
the stress in the body are found in the absence of the flaw, by any suitable method, the stress 
due to the chosen kind of strain nucleus are found, in the same geometry, any number of strain 
nuclei may be distributed through the body without violating the surface boundary conditions. 
The problem therefore remains only to distribute the strain nuclei in such a way as to ensure that 
the crack faces remain traction free, and in order to do this an Integral Equation is established 
which typically must be solved by an appropriate numerical quadrature, this constitute the third 
part of the problem 

In two dimension, a dislocation may be created by first making a slit in an solid, the trajectory 
of slit may take any form, and the dislocation is formed by imposing a constant displacement 
between adjacent points located either side of slit, by inserting or removing the material as nec- 
essary, and then welding the material back together, the constant relative displacement imposed 
is known as Burgers Vector. 

The jargon ’climb’ and ’glide’ dislocation is defined as follows: If the Burgers Vector is 
Perpendicular to the path cut, so that a strip is inserted as shown in figure (1.4), the dislocation 
is said to be ’climb’. If the Burgers Vector is parallel to the path cut, so that there are shear 
displacement, the dislocation is said to be ‘glide’, and the path cut defines the ‘glide plane’. 

1.2 Fundamentals of Distributed Dislocation Technique 

In order to introduce the technique in the simplest possible way consider a plane problem. fig( 1 5) 
shows a plane crack opened by a tensile field. From Buckner’s theorem [ref-12] the solution to 
this problem can be obtained by a superposition of the problems as shown in the figure. These 
are the stresses arising in the uncracked body, and the stresses induced in the unloaded body, due 
to application of equal and opposite traction to those present along the line of the crack in the 
problem. The strategy adopted to generate the corrective traction as shown in the figure is to 
make a fine slit along the line of the crack; the two sides of the cut are then separated by inserting 
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(a) 



Figure 1.4: The creation of an edge dislocation, (a) by an opening displacement (climb) or (b) by 
a tangential displacement (glide) 


material to fill ‘the crack’ as shown, the interior of the real open crack is, of course, empty; the 
inserted material is simply a mathematical device-a means of generating the corrective traction, 
and at the same time simulating separation of the crack-faces. 

It is helpful to visualize the inserted material as a combination of infinitesimally thin strips, 
the first begining at one crack-end and extending to infinity, or a remote boundary, as shown 

in the figure'(l-j3)i; more strips are added as shown in the figure and, finally, strips of material 

■* 

are taken away so that the geometries of figure (1.6); are recovered. The stress field due to the 
complete array of strips may be obtained by a simple addition, or integration, of the solutions 
for single strips. 

The single inserted thin strip of material is the familiar edge dislocation, so that the strain nucleus 
appropriate to plane crack problems, and which we will employ throughout this report, is the 
straight dislocation. 

1.3 Plan of the Thesis 

The second chapter of this thesis covers the Model, Mode-II, problems of a straight and an 
arc shaped crack. Results obtained with the Distributed Dislocation Technique and the results 


6 



0 


0 



Figure 1.5- Application of Buckner’s principle-(a) overall problem (b)stress in perfect body (c) 
separated material 


available in the literature are given in this chapter. This chapter actually proves the validity of 
Distributed Dislocation Technique. Definitions of some basic terms e.g. SIF, Dislocation Density 
etc. are also given in this chapter. Method of handling a singular integral equation is also 
discussed in this chapter. 

Third chapter of this report covers some more complex cases. This chapter covers the problem of 
bi-periodic cracks in a body.effect of periodicity on the SIF is presented with the help of graphs. 
Last part of this chapter covers a circular inclusion problem subjected to uniform pressure. 
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(c) (d) 


Figure 1.6: The Insertion of material between the crack faces by adding and taking away thin 
strips 
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Chapter 2 

Dislocation Techniques 


2.1 Stress Fields 

The stress field of a dislocation located at the origin in a material whose elastic properties are /j 
and v, where b is the burgers vector of the dislocation. 

2.1.1 Stress field of a dislocation with burgers vector in x-direction: 



(a) 



(b) 


Figure 2.1: (a), (b) are Dislocations with Burgers vector in x and y direction respectively 


-fib y(3x 2 + y 2 ) 

27r(l - u) ( x 2 + y 2 ) 2 

(2.1) 

lib y(x 2 - y 2 ) 

2tt(1 — v ) ( x 2 + y 2 ) 2 

(2.2) 
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(2.3) 


lib x(x 2 — y 2 ) 

° xv ~ 2tt(1 - v) \x 2 + y 2 ) 2 

2.1.2 Stress field of dislocation with burgers vector in y-direction: 


fib x(x 2 — y 2 ) 

(2.4) 

2tt(1 - v) (x 2 + y 2) 2 

lib x(3 y 2 + x 2 ) 

2tt( 1 — i/) (x 2 + i/ 2 ) 2 

(2.5) 

lib y{x 2 - y 2 ) 

(2.6) 

27t(1 — u) (x 2 + ?/ 2 ) 2 


2.1.3 Stress Field because of a unit dislocation lying on the boundary 
of a circular inclusion 



Figure 2.2: A dislocation on the boundary of a circular inclusion 


<K 


C 

IT 


/3siny 
2 r 


+ 


smycosy 
2r(l — cosy) 


(2.7) 
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G rB = 7T- [P - (1 - /3)C0S7] 


t C 2 -0(1 + a) , 

^ = 2r7r l + Q -2/3 + (1 “ ^>“ ST 


t _ C fisinj sin'ycos'y 
r 6 tt _ 2r 2r(l — C0S7) 


where, 


2pi(l + a) 


2^(1 - a) 


( 2 . 10 ) 


(2 11 ) 


(*7 + 1X1-^) («a + l)(l-^) v ' 

p 1, r'l and p 2 > ^2 are modulus of rigidity and Poisson’s ratio of the material ‘1’ and ‘2’. a 
and P are the Dundurs parameters given by 

_ /^2(^l + 1) — /ij(/C 2 + 1) 

^2(^1 + 1) + /il(^2 + 1) 

n _ /^2(«l - 1) - Pl(«2 - 1) 

^ 2(^1 + 1) + Pi(«d + 1) 
k is the Kolosov’s const ant :k = 3 — 4^ 

ssuperscript ‘p’ and ‘t’ indicates the stress field because of dislocation with burgers vector in 
radial and tangential direction respectively. 


2.2 Dislocation Density 

The Dislocation Density physically represents the negative of the slope at any point between the 
crack faces. Thus if g(x) represents the amount of opening from the end of a crack then, 


g (x) = - f 
J -a 


( 2 . 12 ) 


dg( 0 


(2.13) 


where ?/>(£) represents the dislocation density at the point £. Note that, since the opening 
displacement is parabolic in form near the crack tips, the density must tend to infinity in a 
Square root singular manner as the crack tip is approached. 
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2.3 Stress Intensity Factor (SIF): 


The Stress intensity factor denoted by K is a measure of the stress singularity at a crack tip 
This stress singularity is square root dependent with the distance from the crack tip. if the 
magnitude of K at a tip is known for a crack which is subjected to a given stress field then K 
remains a measure of how close the crack is to a fracture condition. 

SIF can be computed from the dislocation function density by the following formula' 

K n + ^K I = — - lim - x) [4> x (x) + i^ y {x)\ (2.14) 

— V) x->a 

m terms of non dimensional parameter SIF is given as 

K n + 1 K 1 = — - — r lim 0tt(1 - 77) [ip x (r]) + I'lpyi'q)} (2-15) j 

— V) TJ— »■! 

where Ki\ Stress Intensity factor for Mode-I condition 
I(n ■ Stress Intensity factor for Mode-II condition 

ip x and 'ipy are the dislocation densities for x and y dislocations respectively. I 

I 

j 

2.4 Singular Integral Equation and Gauss-Chebyshev Quadra- j 
ture 1 


Consider the following singular integral : 


S(x) 


1 

7T 



ll)(t)dt 

t — X 


— 1 < X < 1 


(2,16) 


where 

ip{t) = R(t)F(t ), and R(t) = 

let us assume that F(t ) can be approximated to a sufficient degree of accuracy by the following 
truncated series : 

F{t) « Y%.. 0 a l T l {t) 


so now the integral equation becomes 


5(*) 


vn n 1 f +1 Ttfydt __^ n ( ^ 

1=0 v u u - x)^¥ ~ 1=1 1 l ~ l( } 


where the following relation is used : 

1 f +l T t {t)dt _/ Ui-i{x), for i > 0 

n J —i (t - x )VT zr R _ l 0 for * = 0 


(2.17) 
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now foi x — x k 


but for 0 < i < n 


S(x k ) ~ E 1 n =1 a i U l _i(a: J fc) 


= S“ =1 


n(tj - z fc ) 


(2 18) 


and T n (t 3 ) = 0 , C/ n _ x (x k ) = 0 

where T n (x) and £/ n (:r) are chebyshev polynomials defined by 

T n (x) = cos(n6), U n (x) = s -?&+m, cos(9) = x (2.19) 

so our singular integral equation now can be represented in the form of a sum computed at Xk 
points and is given by 

= (2 ' 20) 

where 

T n (tj) = 0, t 3 = cos{-ir(2j - l)/2n), j = 1, ... n / 2i21 \ 

U n -i(x k ) = 0, Xk = cos(kir/n), k — l,...n-l 


2.5 The Mode-I crack : 


2.5.1 Solution using Dislocations: 

Model crack is shown in the fig 2.5.1. The governing integral equation is written below where 
we are computing the stress field at the point x because of a dislocation with unit burgers vector 
situated at the point x'. 


0. 


yy 


/ -j-a 

a y yy (x,xi)'ip y (xi)dx l 

■a 


+ t7°° =0 




To, 


o±(x,x/)i/> x (x/)dxt + / al y (x,xi)ip y {xr)dxr = 0 


rewriting the equation, 
cr*(x,x/) 


^ M x,) = Thfe 


9{ - X,) dx' + r cyyjx, xi) -JM== =dx' + a°°= 0 




X / A 


/■+“ /(*/) 




a / (° 2 - a : 2 ) 


dx 1 — 0 


( 2 . 22 ) 

(2.23) 

(2.24) 


(2.25) 

(2.26) 
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p 



Figure 2.3: Mode-I crack subjected to far field P 


writing the above equation in non dimensional terms 


" +1 . f +l 


-i y/{l~V 2 ) 


g(v)dv+ / 


7-1 \/(l - rf) 


f^drj + a 00 = 0 


,+i a x 


,(S>v) 


f+i n y 


I ’""" g(r>)dn+ I 4&4/M = o 


J-i ytw) ' u dwi 

expressing g(rj) and f{rj) in terms of unknown coefficient and known polynomials, 


g(v) = . f(v) = ^=A T i 


2-1 


sr=i 


’ +1 Vyyi^V) 


1-1 




a-i + 


H ' 1 




^=1 


f,+1 cr^J^ri) 

-1 v 1 — ^ 


i x/ 1-7 ! 

a. + sr=i 1 


L + a°° = 0 


-i x/1 “ 7 2 


71-1(77) 


6 , =0 


applying the quadrature rule 


Z:=i~Xj=i i(ifc)] a, + &. + ff°° = 0 


(2.27) 

(2.28) 

(2.29) 

(2.30) 

(2.31) 

(2.32) 
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(2.33) 


£2=1-21 


n - , 3 =d cr ^A^Vj)T l -i(r] J )] a , + fr, = 0 

applying the no net dislocation condition 


/_ + a a ^ X {xi)dxi = 0 , /_ + ; iJ> y (x/)dTi = 0 


01 


0 (a 2 — XI 2 ) 

in terms of non dimensional parameter 

, '+ 1 p(t?) 






^1=1 




a.=0, S” =1 


f+l T t l ( V ) , 

J " 1 a/W ^ 


6,-0 


s?=.5 [E?.ir, 0, = 0 , E“ =1 J [Sj_ 1 T l _i(^ J )] 6. = 0 

w herel; k — coskn/n , k = 1 , n - 1 

rjj = cos(tt(2j - l)/2n), j = 1, n 

solving equation (2. 31), (2. 32) along with (2.37) for the unknowns 

K i = v /2?r ( 1 - *7) [A ( 7 )] 


70 = 



r lim ->/27r(l - ??) 

I/) 7J-+1 


SkMlzlM 

%/! “7 2 


70 


2a(l — i/) 



2.5.2 Analytical result: 


77 / = Py/wa 


(2.34) 


(2.35) 


(2.36) 


(2.37) 

(2.38) 

(2.39) 


(2.40) 


2.5.3 Comparison between Computational and Analytical Results 

For (j, = 80 Gpa, ^ = 0.3, P = 80 MPa, 2a = 20 mm 
7f / [computational) ~~ 14.1796308 MPa \[rfi 
1 (analytical) = 14,1796308 MPa 

while Chan et al. (1970) has obtained the results by FEM with an error of about 11% with 
coarse mesh, and about 5% with fine mesh. 
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2.6 The Mode-II crack : 

2.6.1 Solution using Dislocations: 

Modc-II cxack is shown in the figure. The governing integral equation for this mode is same as 
that for model 


P 



Figure 2.4: Mode-II crack subjected to far field P 


/ ■+a r+a 

a* y (x, xr)ip x (xi)dx' + / a” y (x, x/)ip y (x/)dx' = 0 

-a J -a. 

/ -ha P+ a 

cr xy {x, xi)ijj x (xi)dxi + / a y xy (x,xf)ip y {x/)dxf + r°° = 0 
J ~ a 


AM = = 7 

y (a 2 - a; 2 ) V (a 2 - 


rewriting the equation, 



(2.41) 

(2.42) 

(2.43) 


(2.44) 
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a x (x,xr)- 


(a 2 — x 2 


/‘+a 

--dx'+ o y xy {x,xi)- 

J —a 


wuLing the above equation in non dimensional terms 


V 2 -^) 


dx' + r°° = 0 


(2 45) 




representing g(g), f(t]) in terms of unknown coefficient and known functions 


--9(v)dr]+ [ 

) J - 1 


+1 a xy^,V) 


f{v) + r °° = o 


(2.46) 


(2.47) 


g(v) = Z: =1 aX-i(v),f(v) = 


(2.48) 


writing the above equation in non dimensional terms 


a yy(^vW'n)dv+ / a y ((,g)^(ri)dg = 0 


+ r°° =0 


/* ^x V (^v)^{v)dr] + J *%y{Z>ri)1>{y) 


applying the quadrature rule 


(2.49) 


(2.50) 


2,"=i [<*(&. %)r.- t fe)] 0, + L ^ E” =1 [a*,(&, n,)T,-i(%)] b, = 0 

77 / 1 1 


(2.51) 


O. + SJ=^EJ.iK(&. ■?,)?.., (%)] 5. +T» = 0 (2.52) 


applying the no net dislocation condition 


o-|-a r+a 

i/) x (xi)dxt = 0, / ip y (xt)dxi - 0 


(2.53) 


J —a ^/(a 2 -x/ 2 ) 
in terms of non dimensional parameter 


da;/ = 0, 


nda;/ = 0 


+i _iE=d,=o, r l -m f *,=o 

i 7-1 A/ 1 -/ 7 2 


(2.54) 


(2.55) 


£?=, [ r °i= o ' s -' 


Tt-i(rj) 


drj\b z = 0 


(2.56) 
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(2.57) 


^ I=l 77. SJL 1 — [52” =1 T)_ i ( 77 j )] b, ~ 0 

I V 


wheic £a = coskn/n , k=l,n-l 
rjj - cos(n(2j - i)/2 n) j=l,n 

solving (2 o2), (2.o3) along with (2.58) for the unknowns, and computing the stress intensity 
factor with flic help of formula (2.61) given below, 


K, 


n — 




(2.58) 


K u = 




2a(l - v) r£i 


lim 1 / 271 - (1 - rj) 


T^ =l aX-\{r)) 

V 1 - 7 ) 2 


Kir — 




2a(l — 1 /) 




2.6.2 Analytical Result: 


K 11 = P 1/702 


(2.59) 

(2.60) 


(2.61) 


2.6.3 Comparison between Analytical and Computational Results 

For /i — 80 Gpa, // = 0.3, P = 80 MPa, 2a = 20 mm 
K a (computational) — 14.1796308 MPa \frn 
Kn(analyttral) = 14.1796308 MPa y/m 


2.7 Arc shaped crack : 


2.7.1 Solution using Dislocation: 


The governing integral equation for a arc shaped crack is written below.where we are computing 
the stress field at a point 9 because of a dislocation located at an angle <j>, and the crack is 
subjected to uniform pressure. 






<%o(PMM) d 4> + 



0*r(0, + °Vr(<£) = 0 



aj, 0 (9, 0)0/0)# + a r e{<t>) = 0 


where superscript ‘p’ and ‘t’ are used for the dislocations with Burgers vector 
tangential direction respectively. 

2a:Included arc angle. 



9 {4) 

\/(Q! 2 - 0 2 )’ 


0 t(0) = 


f{4>) 


in 


(2.62) 

(2.63) 
radial and 


(2.64) 
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Figure 2.5: Arc shaped crack subjected to uniform pressure 


rewriting the equations 


<( 0 ' ,l, '>-7=rh? d4 ' + 


/ +a 

■a 


/ a 2 - (j > 2 


d(j) + eYr(^) = 0 




4=~d<j) + g t0 {<1)) - V 

J-a A “ 0 


in terms of non dimensional parameters 


fr(£> r /) 


y 

I rto&v)- 


'1 _ ? ^2 


=d?7 + cr rr (£) = 0 


rtZty + ffrfl^) = 0 


(2.65) 


( 2 . 66 ) 


(2.67) 


(2.68) 


jy rlM -0L ^ + £ (2 ' 68) 

9 („), /(„) can be expanded in terms of chebyshev polynomials with unknown eoefficients a, and 


g (n) = ^yo.T^Jiv) = 


(2.69) 


■ +1 i&MlU (,) J O. + Ef-! / + ‘ #= r <- fo)*> 6 ‘ + ' *- = 0 (2 - 70) 

-1 a/1 - v 2 n _1 V / 
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S T! 

7=1 


+ 1 


<e(^) 

x/T 


-i \/i -r}‘ 
applying the quadrature rule, 


T t -i(ri)dri 


a i "h sr = i 


r< 
j - 1 ^ 


<T t rO& r l) 

V 1 -v 2 


T l - l {r])d'i] 


b t + cr r o — 0 


s;=i -s j ; =i K r (a, o,j + £r =1 -£? =1 KV(a-, ^-iM k + a„{^) = o 

/t' Ti 

s;u !!j)r,_i(rjj)] a, + E» = ds» , ^T, -,(»,)] b, + a T „(( t ) = 0 


n 


no net dislocation condition in global (x,y) set 

r>4-Q; 




il>x(4>)d<l> = o 


applying the transformation, 


(2 71) 

(2.72) 
(2 73) 

(2.74) 


Y 



X 

Y 


sin0 cos0 


r 

cos0 -sin0 


t 


Figure 2.6: Rotation Matrix 


or 



[^ p ((l>)s%n{(f>) + ^i(tf>)cos(^)] # = 0 




9 ^ - sin(i j)) + — 7====== cos ( 4 >) d(f> = 0 


l-a 


vV^ 2 ) 


(2.75) 

(2.76) 
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in tonus of non dimensional parameter 


+1 

i 


3 ^ m(na) + — ~ L=cos(ria) 




/m2 


yfi - ?] 2 


dr] = 0 


V" 
^i- l 


11 r /;-i (»/) , ' 

—===tiin(rja)dri 
i v/ 1 ~ '0 2 


«. + 


+i 


Tt-ifo) 


=cos(rja)dr] 


k = 0 




-1 yl~- 7/ 2 
7T 

[^7=1-^*— 1 ( T ?j) s * n ( 7 7j Q£ )] a t + S" =1 — [S”_ 1 Tj_ 1 (?7 j )co5(7] j q;)] b t — 0 
for the no net dislocation in y-direction 

/ +Q 

%{<j))d(j) = 0 


-I- tv 


[i)p{(j))cos{(j>) - ijji((j>)sm((j))] d<$> = 0 


following the similar steps 


^"-i” , 7;_i ( t/j ) cos(r 7 jQ!)] o t - [S^ =1 T,_ 1 (r 7 J )stn(j]i J a)] 6, = 0 

where = coskn/n, k=l,n-l 
Vj = con (it (2 j - l)/2»), j=l,n 
in terms of non dimensional parameter SIF is given as 

K n + iKi = lim /Ml - 77 ) [M/ + M/)] 

/Oi + iKr = 2 Q( i~ ~ 1 / ) ^ F ^ flt + 


(2.77) 

(2.78) 

(2.79) 

(2.80) 
(2.81) 

(2.82) 


(2.83) 

(2.84) 


2.7.2 Analytical Result: 


, cos f 

K Ia — avirrsma 2a 

A l+sin 2 ^ 

(2.85) 

, smf 

Kn, - O'jir Bttj + - sm2f 

(2.86) 


2.7.3 Comparison between Computational and Analytical results 

For 11 = 80 GPa, u = 0.3, a = IMPa, r = 0.05 m, 2a = 180 degrees 
Kj, {computational) = 0.186833041 MPa 
K I(analvtical) = 0.186833041 MPa 

r-;-~;,,r5pr 

SfiTTf^T-w * A . 0.3 62E 
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2 


Figure 2.7: Graph Showing the variation of Ku vs N: 


Graph Showing the variation of SIF vs 2 a 





Chapter 3 

Complex geometries 

3.1 Bi-Periodic cracks: 


c 



a 


Pig ur e 3.1; ]3i periodic cracks in an infinite body with periodicity alpha in x-direction and Beta 
in Y-direction: 


Fig(3,l) shows bi-periodic cracks in an infinite body with periodicity a and /? in X and Y 
direction respectively subjected to far field stress a an r. The stress field of dislocation for 
this case is taken from reference- [7]. Graph shows the variation of stress intensity factor with 
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parameter a and (3 respectively. 


3.1.1 Array of cracks with periodicity Alpha in X-direction: 

F'tr>(3 2) shows the periodic array of cracks with periodicity a. in x-direction (when f3 is much 
gi cater than a ). variation of SIF with respect to different values of a and jB is given in the 
graphs 



Figure 3.2: Array of cracks in an infinite body with periodicity a in x-direction where ft >> a 


3.1.2 Array of cracks with periodicity Beta in Y-direction: 

Fig(3.4) shows the periodic array of cracks with periodicity a in x-direction(when a is much 
greater than (3 ). variation of SIF with respect to different values of a and /3 is given in the 
graphs. 

actual data for KI are taken by Fourier Transform (Fichter) , Expansion of Complex Poten- 
tial(Isida) Asymptotic Approximation (Polynomial of degree ten) (Benthem) . Accuracy of these 
data is about 2%. 

3.2 Circular inclusion with interface 

In composite materials containing fiber or particulate reinforcement the interface separating ma- 
trix from inclusion is widely believed to be a dominant influence affecting the overall stiffness and 
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K 

Graph between /K Q And Alpha 


’beta3 dat’ 

’beta4 dat’ 

’beta8.dat’ 

’betal6.dat’ 
’beta32 dat’ ~ ~~ 
sqrt(x" i 'tan(pi/x)/pi) n 



Figure 3.3: Graphs showing the variation of Kj/Kq vs periodicity in x direction i.e. Alpha for 
different values of beta 


Graph Showing variation of 


vs ot 


Alpha 


’k2b4 dat’ 

’k2b8.dat’ - 

’k2bl6 dat’ 

’k2b32 dat’ ■ 

sqrt(x*tan(pu'x)/pi) a 


Figure 3.4: Graphs showing the variation of K n /K 0 vs periodicity in x direction i.e. Alpha for 
different values of beta 






damage tolerance characteristics of the composite. Damage accumulation in a composite often 
depends on the character of the mechanical response of the interface. Hence accurate assessment 
of the overall stillness characteristics of a composite containing fiber bonded to the matrix is an 
extremely important in the attempt to obtain improved composite performance. 

Standard material models require for their implementation the solution to a solitary mclusion- 
iiit.erface-mal.rix system the problem of obtaining effective material properties for a composite is, 
to a large extent, dependent on the availability of the elastic field solutions to the solitary problem. 


3.2.1 Interface Law 

The interface separating the matrix from inclusion is characterised by a constitute equation 
relating the components of interface traction to components of interface drsplacement jump, 
only normal inter-facial traction develops from the normal component of drsplacement jump ren, 
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p 


Figure 3.(>: Graphs showing the variation of KI/KO vs periodicity in y direction i.e. Beta for 
different values of Alpha 


owing to rotational symmetry, the force law assumes the form 

*5* (n) = /(K])n 

where S( n) is the interface traction vector acting on the interface and f is a scalar function of 
radial interface displacement jump [u r ]. 

3.2.2 Problem Formulation 

consider a circular inelusion-iuterface-matrix system subjected to uniform pressure a , where V 
is the radius of the circular inclusion, 

The governing integral equations for a circular inclusion problem is given below. 

/*2 7T 

/ 0 •';,((?, + / 4r(0> + °rr( 9 ) + a rr (u r (9) , U t (0)) = 0 (3 1) 

Jo Jo 

I ^ + f °r0^' + °reM 9 )M 9 )) = o M 

Jo Jo 

where superscript ‘a’ indicates the stress field because of applied pressure, and c indicates the 

stress field because of inter-facial bonding,*. e. cohesive zones. 

u r (0) is used for jump-in- displacement in radial direction at the point 9 . 
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a 


Figure 3.7: A circular inclusion subjected to uniform far field pressure©- 

similarly u t (6) indicates jump -in- displacement along tangential direction at the point 9. 
writing the above equations in non-dimensional terms and making use of interface law , 

/•+ 1 r + 1 

a J ^ <•(£, v)Mv)dv + a J 1 a rr(^ n )Mv)dv + <(£) + Ku r (0 = 0 

/'+i /*+i 

“ J , + a J 1 = o 

wliorc, = Kur[(), o?«(«r(Oi “iR)) = ; 4“iR) 


(3.3) 

(3.4) 


AT : Interface Modulus 


u x (£) = a 0 + aj il> x (Ti)dri 


where, ao : Jump - in — displacement at (r,o) 

u x (f) = «o + a J bP p{rj)cos9 - i/i ( (?])s*n0] dp 

expressing dislocation densities in terms of unknown coefficients and chebyshev polynomials 

MV) = X%ZZa p T p (r,)Mv) = Eg b p T p (p ) 
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»'.r(0 -- <h\ -!- (>mJa p 

applying tin' quadratino rule, 


J Tp(rj)cosQdj) 


<*££ X 


u-i 


T p (r})sm$dr j 


whci o 


similarly 


vAQ ■- 


<!;) - «o -I- (i L j a p \H ((, rjj)T p (rij)cos9j 

~ a ^p=u a p^]=\H{^ty)T p {r] :) )sind J AjT-^ 

II (f n) - S °> for £ < v 

\ 1, for -1<?7<£ 

%(£) = &o + a J i> y {r))dr] 


(3.5) 


whore, h a : Jump - ;in - displacement at (r,o) 

K 


(£) =io + « j [ip p {ri)szn8 + ip t (7])cos9] dp 


expressing dislocation densities in terms of chebyshev polynomial and applying the quadrature 
rule, 

u, 


<»(0 = bo + o.Z;Z fa, S^m rhKMsmS, 

n,)T r (v,)cos6, 073^ 


(3.6) 


u r (6 = u x (g)cos9t: + u y {i)smB^ 

applying the transformation 

Ur(0 ~ (IqconO + bosinO -I- u S()_ q — S^_ ^ id ( ^ , Pj)Tp{Vj) yjl — r}jCos(9^ — f?j) 


(3.7) 


+a^Z 2 0 b p -^M>r} J )TMJ T-rgMOt ~ * 3 ) 

Tl 

ut{0 = -«nfl € ii*(f) + cosOtUyiCj 
substituting the values of u w and % in the expression of u t 

u^) ~ -ao sin0( + b 0 cosd{ ~ aEpZl a r ~ (£ > ^ ) (% ) ~ %) 


(3.8) 


(3.9) 


(£, %)T P (%) ^Jl-r^cos ^ ~ fy) (3.10) 
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cos 9 sin0 u x 


sin0 cos0 u 


Figure 3.8: Rotation Matrix 

substituting tin* values of u r and u t in the governing integral equations and making use of 
quadrature rule, after readjusting the terms the equations become 

>? K(& %) + nWf - <t)l 

I M” ffc,,” >™ i d 1 - t)j« K(6. 'h) + M(&. rjjjsmlfl; - »,)) 

-VKcosO^ + Ksin9^b 0 = -a“ r (^) (3.11) 


V»» ?t U« 7> 

* (I 'gj * M !> v/j . 


^ - KH{Z k ,ri 3 )sin(9z - 0,)] 


i >;;; y./;, (?/.,•) K(&, »&•) + - 0>)] 

-Ksind^ao -4- I<cos9^b 0 = -a“ g (£k) 

imposing the no net dislocation condition in global (x,y) coordinate system 

/ Mv)dv = 0 


(3.12) 


[cOfltf^fo) - «nfyV»t (»/)] dr ) 



applying the quadrature rule, 


S p=o a p~^i T p(r}j) c os8j yjl- rf - E£ =0 \ - E” =1 T p ) sind 3 yfl-r\ 2 = 0 (3 13) 

7 % 


similarly for y dislocation, 

r + 1 


-+1 


'-i 


i>{ri)dri= / [sin0jil> p (n)+cos6jil>t(v)]dri = O 


applying the quadrature rule, 


s;u 2 ^~^ =1 T p ( % ) 5 m^ V 1 - ^ + = 0 (3 14) 

cquation(3.11), (3.12), (3.13) and (3.14) now solved for unknowns. once the unknowns are deter- 
mined one can compute the jump-in-displacement at any point on the circular inclusion by using 
the equation (3.7)and (3.9)respectively. 


3.2.3 Analytical Result 

For the plain strain problems the displacement vector in radial direction U r is given by: 

U r = ClT+~ (3-15) 

but if the material is homogeneous without any voids or hole then above equation reduces to, 

U r = dr (3.16) 

jump in displacement at the interface is given by the expression, 

u r = lim U r — lim U T (3.17) 

r~ ir+ r-+r - 

to compute the unknowns we make use of the boundary conditions of traction, continuity and 
arc given by 

lim oy = f(u T ) 

r— >r+ 

lim U T = f[u r ) 

r—>r~ 

lim ay = cr® r (3.18) 

r-+oo 

where / is a scalar function of radial interface displacement jump, making use of above boundary 
conditions for the unknowns, the jump-in-displacement is comes out to be, 

^^ = 2(1 - v 2 ) [cr - f {u r )] (3.19) 

r 

where the matrix and inclusion are assumed to have same elastic properties. 
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3.2.4 Comparison between analytical and computational results 


for, E = 18.2 MPa, r = 1mm., a = IMPa, u = 0.3 
u r (computational) = 0.090909mm 
u r (analytical) = 0.090909mm 



Chapter 4 

Conclusion and Future Scope 


One of the advantages of this technique is that it is possible to determine explicitly the state of 
stress induced by a dislocation in a range of different geometries in an infinite plane, a half plane, 
a strip or circular inclusion etc. it is possible to distribute any number of dislocations within 
the body in question, whilst preserving the boundary condition on the surface of the object, the 
problem therefore reduces to one of determining the distribution of dislocation which generates 
tractions, equal and opposite to those produced by external loading. 

In chapter-2 the SIF computed for Model, and Mode-II cases for a straight and a circular crack 
matches exactly with the available results given under Analytical Section. In case of bi-periodic 
cracks the computed SIF matches exactly with the available results.in case of bi-periodic cracks 
the the SIF is found to be exactly Psfrfd for large values of parameter a and /?, which is the case 
of a single crack present in an infinite body. 

The inclusion problem also gives satisfactory resuits.The parameter ‘K’ in the inclusion problem 
determines the kind of geometry being considered. for K = 0 we are characterizing a hole in the 
body. The opening computed for this case also matches with the opening obtained with the help 

of rcf-[l]. 

From the previous discussion it is clear that Distributed Dislocation Technique can be used to 
find solution related to fracture mechanics with full accuracy. 


Its future scope are, 


• Using Distributed Dislocation method one can solve other crack problems involving com- 
plicated pattern of cracks in the body. e.g. Staggered cracks in a body. 

• Inclusion problem can be extended to other type of geometry probiems.one can use this 
method to find the solution for elliptical or other type of arbitrary shaped inclusion. 
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Inclusion pioblem can be used to find the effective properties of a material containing voids. 


a 


• Inclusion problem can be solved for the case where matrix and inclusion have different elastic 
properties. 

• Non linear cases, where the scalar function of radial displacement jump is not linear depen- 
dent on the radial displacement jump, can also be solved. 
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